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Recent  theoretical  calculations  of  electron  density  profiles  within  the  high  latitude/polar  cap 
ionosphere  demonstrate  that  the  diurnal  foF2  variation  observed  at  Thule,  is  controlled  by  the 

plasma  convection  pattern  and  the  associated  drift  velocities.  The  model  calculations  for  Bz  <  0 

and  Bz  -  0  show  factor  2  to  3  differences  in  ftnax  over  Thule,  supporting  the  stated  importance  of 

convection  pattern  and  velocity  measurements  for  the  modelling  of  the  high  latitude  ionosphere. 
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ABSTRACT 

Plasma  convection  measurements,  using  Oigisonde  ionospheric  sounders,  have  been  conducted  in 
the  central  polar  cap  at  Thule,  Greenland  (86°  CGL)  and  more  recently  at  Qaanaaq,  Greenland  (87° 

CGL),  in  the  oval  at  the  AFGL  Goose  Bay  Ionospheric  Observatory  (65°  CGL),  and  at  suboval  latitudes 
at  Argentia  NAS  (57°  CGL).  The  plasma  convection  or  ionospheric  drift  measurements  conducted  at 
Thule  and  Qaanaaq  during  campaigns  from  Winter  1983/84  to  present  provide  evidence,  that  antisunward 
convection  dominates  in  the  polar  cap  with  velocities  typically  between  300  and  900  msec*1.  Drift 
shoals  were  observed  during  periods  of  arc-transition  (quiet  magnetic  conditions).  Observations  of 
the  plasma  drift  at  Goose  Bay  show,  as  expected,  a  drift  reversal  from  westward  to  eastward  around 
midnight  COLT,  indicating  the  transit  of  Goose  Bay  from  the  evening  to  the  morning  convection  cell. 
Observations  at  Argentia,  typically  a  suboval /trough  station,  provide  evidence  under  magnetically 
disturbed  conditions  for  the  midnight  reversal  of  the  antisunward  flow  pattern.  However,  the  data 
are  less  consistent  under  magnetically  quiet  conditions.  The  proximity  of  the  station  to  the 
boundary  between  corotating  and  convecting  plasma  may  at  times  affect  the  consistency  of  the 
measurements. 

Recent  theoretical  calculations  of  electron  density  profiles  within  the  high  1 atitude/pol ar  cap 
ionosphere  demonstrate  that  the  diurnal  foF2  variation  observed  at  Thule,  is  controlled  by  the 
plasma  convection  pattern  and  the  associated  drift  velocities.  The  model  calculations  for  Bz  <  0 
and  Bz  «  0  show  factor  2  to  3  differences  in  Nnax  over  Thule,  supporting  the  stated  importance  of 
convection  pattern  and  velocity  measurements  for  the  modelling  of  the  high  latitude  ionosphere. 

INTRODUCTION 

The  I  IT  dependence  of  the  winter  polar  cap  ionosphere  (Sato  and  Rourke,  1964,  Buchau  et  al., 

1985),  the  longitudinal  variation  in  F-region  densities  in  the  nightside  auroral  zone  (de  la 
Beaujardiere  et  al.,  1935)  and  the  formation  of  the  main  F-layer  trough  are  either  entirely  or  to  a 
large  extent  the  result  of  large  scale  plasma  convection  across  the  polar  cap.  (Knudsen,  1974,  Knudsen 
et  al.,  1977,  Watkins,  1978,  Sojka  et  al.,  1981)  and  of  the  displacement  of  the  geomagnetic  and 
geographic  poles  (Sojka  et  al.,  1979). 

The  existence  of  polar  plasma  convection  has  been  inferred  from  electric  field  measurements  from 
satellites  (Caufman  and  Gurnett,  1982),  or  by  directly  measuring  the  convection  with  satellite-borne 
plasma  drift  measurements  (Hanson  and  Heelis,  1975).  Ground  based  high  latitude  plasma  convection 
measurements  have  been  conducted  during  the  decade-long  measuring  campaign  of  the  Chatenika  Incoherent 
Scatter  Radar  -  ISR  (Foster,  1983)  and  looking  northward  by  the  Millstone  Hill  I  SR  (Evans  et  al.,  1980). 
The  large  body  of  plasma  convection  data  shows  that  in  its  simplest  form  the  convection  occurs  in  a 
two  cell  pattern.  Models  of  this  two  cell  convection  pattern  developed  by  Volland  (1975)  and  improved 
by  Heelis  et  al .  (1932)  have  been  used  extensively  to  model  the  ionosphere  under  the  influence  of  the 
convection  (above  references).  Figure  1  from  Heelis  et  al  .  (1982)  shows  the  two  cell  configuration, 
with  antisunward  p'asma  flow  across  the  polar  cap,  and  sunward  return  flow  along  the  morning  and 
evening  flanks  of  the  auroral  oval. 
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Figure  1.  Polar  Plasma  Convection  Pattern  (Heelis  et  al.,  1982)  in  CGL/CGLT  coordinate  system. 
The  locations  of  the  drift  measuring  stations  Qaanaaq,  Goose  Bay  and  Argentia  are  indicated  by 
dashed  circles. 


Past  observations  and  modelling  efforts  have  shown  that  knowledge  of  the  convection  pattern  is  a 
prerequisite  for  accurately  modelling  the  ionosphere  over  the  full  high  latitude  region.  The  two 
cell  p actern  occurs  most  clearly  when  the  Interpl anetary  Magnetic  Field  (IMF)  has  a  southward  component 
(B2  <  0),  while  the  exact  flow  geometry  is  controlled  by  By  (Heelis  1984).  For  Bz  >  0,  dusk  to  dawn 
polar  cap  electric  fields  have  been  measured,  which  cause  sunward  plasma  flow  across  the  central 
polar  cap  (Burke  et  al.,  1979)  and  leads  to  an  interpretation  of  the  existence  of  a  four  cell  con¬ 
vection  pattern.  Satellite  magnetometer  measurements  (Zanetti  et  a).,  1984)  suggest  a  three  or  four 
cell  pattern  depending  of  the  sign  of  By.  Lassen  (1979)  interpreted  the  polar  cap  auroral  arcs, 
observed  during  magnetically  quiet  periods,  as  being  the  result  of  a  shear  between  a  third  central 
polar  cap  cell  occuring  either  at  the  edge  of  the  morning  (when  By  <  0)  or  the  evening  cell  (when 
By  >  0).  Recent  analysis  of  DE2  data  (Heppner  and  Maynard,  1987)  indicates,  that  it  is  possible  to 
explain  the  data  for  the  full  range  of  Bz  and  By  with  a  basic  two  cell  pattern.  For  the  Bz  >  0 
condition  this  approach  results  in  substantial  distortion  of  the  two  cells,  in  order  to  explain 
sunward  flow  in  the  central  polar  cap,  but  the  resulting  pattern  agrees  with  observations  in  other 
regions  of  the  high  latitude  ionosphere. 

In  this  paper  we  describe  a  ground  based  coherent  HF  radar  (digital  ionosonde)  technique,  which 
allows  continuous,  relatively  low  cost  measurement  of  the  plasm  convection  at  selected  sites.  Early 
observations  (Reinisch  et  al.,  1987)  show  good  agreement  between  the  measurements  and  the  major 
features  of  the  plasma  convection.  The  analysis  presented  here  indicates,  that  not  only  the  major 
features,  but  also  details  such  as  shears  and  drift  reversals  in  response  to  quiet  magnetic  conditions 
are  observed.  The  measurements  suggest  that  the  continuous  monitoring  of  the  convection  pattern  with 
a  ground  based  technique  will  be  possible. 

SPACED  ANTENNA  DOPPLER  DRIFT  TECHNIQUE 

The  drift  technique,  which  has  been  developed  for  a  digital  ionosonde,  the  Digisonde  (Bib!  and 
Reinisch,  1978),  is  based  on  a  technique  using  spectral  analysis  of  ionospheric  drift  data,  developed 
by  ofister  (1974  a,  b).  This  techni que  has  been  refined  for  automatic  data  processing  (Dozois, 

1983),  and  is  an  integrated  capability  of  all  modern  Digisondes  (Reinisch  1986). 

The  drift  measurements  are  conducted  using  phase  coherent  pulse  transmission  on  two  HF  fre¬ 
quencies,  which  are  selected  to  sample  desired  heights  in  the  F-region.  These  vertically  transmitted 
HF  waves  illuminate  a  large  area  of  several  hundred  kilometers  diameter  in  the  F-region;  an  array  of 
4  or  7  antennas  (depending  on  the  site)  deployed  in  a  100  m  baseline  triangle  receives  the  signals 
reflected  from  the  ionosphere  (the  technique  is  equally  applicable  to  E-region  studies).  The  indi¬ 
vidual  antennas  of  the  receiving  arrary  are  multiplexed  at  the  pulse  repetition  rate  (200  Hz).  The 
time  series  received  at  each  of  the  antennas  is  Fourier  transformed  in  real  time,  resulting  in  four 
(or  seven)  complex  spectra  at  the  end  of  each  measurement  period. 
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Cross-correlation  P  of  complex  Doppler  spectrum  F  from  4  antennas: 

|  -  spectra!  kne 

P  -££  F.  F*.  exp -ih[a -S;.] 

k.8  j  •»  Jl  )I  J  J  a  -  antenna  posit  ton 

yields  angle  of  arrival  of  each  reflected  spectral  line. 


Least  Squares  Fit  for  Drift  Velocity  V: 


a  |  *s<Y-S34  tc>’ 

£-  5*. 


a.-  unit  position  vector  ol  reflector  • 
a  Doppler  frequency  of  reflector  s 
f  -  sounding  frequency 
weighting  factor 
c  -  speed  of  light  in  vacuo 


Figure  3.  Drift  Skymap  showing  12  sources  with  positive  Doppler  and  5  sources  with  negative  Doppler. 

As  will  be  shown  later,  the  size  (or  diameter)  of  the  convection  pattern  and  the  actual  convec¬ 
tion  velocities  are  important  quantities  in  modeling  the  polar  ionosphere.  Three  Digisonde  equipped 
stations,  namely  Qaanaaq,  Greenland  (87°  CGL),  Goose  Bay,  Labrador  (65°  CGL)  and  Aryentia,  NAS, 
Newfoundland  (57°  CGL)  provide  a  first  step  at  defining  these  parameters,  at  least  for  a  large  part 
of  tne  day.  Their  respective  locations  over  a  full  day  are  indicated  in  Figure  1,  by  dashed  circles 
marked  Q  (Qaanaaq),  G  (Goose  Bay)  and  A  (Argentia).  The  figure  shows  that  for  the  conditions  repre¬ 
sented  by  this  pattern,  Qaanaaq  is  fully  embedded  in  the  anti-sunward  flow  in  the  central  polar  cap, 
while  Goose  Bay  is  under  the  influence  of  westward  (sunward)  flow  before  midnight,  and  observes 
eastward  (sunward)  flow  after  midnight.  In  the  example  chosen,  Argentia  would  be  outside  the  high 
latitude  convection  and  would  observe  mid-latitude  co-rotating  plasma.  Under  more  active  conditions 
(in  response  to  8Z  <<  0),  the  convection  pattern  expands  in  diameter,  bringing  Argentia  under  its 
influence  for  part  of  the  night,  and  extending  the  duration,  during  which  Goose  Bay  observes  sunward 
return  flow.  Figure  4  shows  in  a  CG  direction/CG  local  time  diagram  the  approximate  drift  direction 
changes,  which  over  a  24  hour  day  would  be  observed  at  Qaanaaq  (top)  and  Goose  Bay  (bottom),  assuming 
the  convection  pattern  shown  in  Figure  1.  The  top  graph  shows  the  360°  rotation  in  24  hours  of  the 
anti-sunward  drift  at  Qaanaaq,  while  the  bottom  graph  shows  the  westward  drift  before  midnight  and 
the  transition  to  eastward  flow  through  south  at/or  near  midnight  CGLT  observed  at  Goose  Bay.  The 
coordinate  system  of  Figure  4  has  been  chosen  for  the  following  discussion  of  the  drift  observations. 
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Figure  4.  Expected  diurnal  change 
of  the  drift  direction  (assuming  the 
convection  shown  in  Figure  1)  for 
Qaanaaq  and  Goose  Bay  in  a  CG  direction 
and  CG  local  time  coordinate  system. 
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DRIFT  MEASUREMENTS 


a.  Tkylp  -  Magnetically  Active  Period,  1983/1984  Measurements 

Several  winter  measurements  campaigns  were  conducted  by  the  Airborne  Ionospheric  Observatory, 
while  on  the  ground  at  Thule  AR,  Greenland  (86°  CGL)  to  study  the  structure  and  dynamics  of  the  polar 
cat'  iooosrhere.  During  these  measurements  the  antisunward  convection  of  large  non- oral  plasma 
patch“s  was  observed  with  an  All  Sky  Imaging  Photometer  -  ASIP,  (Weber  et  al.,  1984),  and  the  UT 
dep°"‘->ncn  of  the  maximum  electron  densities  in  these  patches  was  documented  (Ruchau  et  al.,  1983). 
For  t"e  campaiqns  in  December  1983,  January/February  1984  and  March  1985  a  four-antenna  array  was 
deployed  nnar  the  aircraft,  to  provide  measurements  of  the  plasma  drift  independent  of  the  optical 
observations  (the  patches  observable  by  the  ASIP  are  generally  absent  in  the  00-12  UT  time  window). 

While  the  spaced  antenna  Doppler  drift  measurements  were  developed  at  the  AFRL  Goose  Ray 
Innosp"eri c  Observatory,  the  best  evidence  that  this  technique  actually  allows  the  measurement  of 
nlasr"‘.  drifts  was  obtained  from  the  initial  central  polar  cap  observations  (winter  83/84).  Figure  5 
shows  the  drift  results  in  the  data  format  discussed  in  Figure  4,  for  measurements  taken  during  four 
disturbed  to  moderately  disturbed  days  between  30  January  and  4  February  1984  (£  Kp  for  30, 

31  Tan.,  1,  2,  3,  4  Feb  was  34,  29-,  23+,  22-,  23,  41+  respectively).  The  data  show  predomi nantel y 
an+i-su"ward  drift  (the  anti-sunward  direction,  rotating  through  360°  in  24  hours  is  indicated), 
corroborating  coincident  optical  observation  for  3/4  February  1984  (Weber  et  al.,  1986),  and  earlier 
reports  derived  from  optical  data  of  antisunward  plasma  transport  under  disturbed  conditions  (Weber 
et .  al.,  1984).  The  observed  velocities  shown  in  the  top  panel  of  Figure  5  (median  velocity  and 
range1  fluctuate  between  150  and  900  msec”',  well  within  range  of  previous  satellite  measurements  and 
in  general  agreement  with  the  optical  observations .  With  the  confidence  in  the  technique  gained  from 
t'-’nse  ■"‘•'’asurnmpnts ,  permanent  digital  sounding  capabilities  were  deployed  at  the  Qaanaaq  Geophysical 
Oh s° r v a t o ry  of  the  Danish  Meteorological  Institute,  in  Northern  Greenland  (87°  CGL),  and  at  a  suboval 
station  .'Argent ia,  Newfoundland,  58°  CGL).  With  the  AFGL  Goose  Bay  Ionospheric  Observatory  (65°  CGL) 
these  station  form  a  meridian  chain  along  approximately  20°  East  CG  Longitude. 
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Figure  5.  F-rngion  drift  direction  (bottom)  and  velocity  (top,  showing  mean  velocity  and  range  of 
observations!  tahen  during  four  moderately  disturbed  days  in  Janurary  and  February  1984  at  Thule  AR, 
Greenland.  The  antisunward  direction  is  indicated  for  reference. 

b.  Thule  .  Magnet. i cal  ly  Ouiet  Period 

The  earlier  polar  cap  campaigns  had  already  established  that  sunaligned  polar  cap  arcs,  espe¬ 
cially  suhvisual  F-layer  arcs,  populate  the  polar  cap  under  magnetically  quiet  conditions  (Weber  and 
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Buchau,  1981,  Buchau  et  al.,  1983).  Reitt  et  al.  (1978)  had  shown  that  current  continuity  at  plasma 
flow  shears  requires  upward  (or  downward)  Birkeland  currents,  therefore  producing  (sub)  visual  signa¬ 
tures  i.e.  auroras  in  the  region  of  upward  currents.  Conversely  Lassen  (1979)  concluded  that  polar 
cap  arcs  should  be  plasma  shear  demarcation  lines.  The  analysis  of  Sondrestrom  1SR  data  in  the 
context  of  ASIP  arc  observations  has  confirmed  this  hypothesis  (Carlson  et  ai.,  1983). 

The  Thule  drift  data  taken  on  9/10  December  1983  cover  a  nine-hour  period  at  the  end  of  a  27 
hour  long  quiet  period  (9  December  Z  Kp  =  5+),  during  which  subvisual  and  visual  arcs  typical  for 
quiet  periods  were  observed,  occupying  major  parts  of  the  sky  and  drifting  from  dawn  to  dusk  from  1900 
to  2300  UT  on  9  December  1983.  The  concurrent  drift  measurements  are  shown  in  Figure  6,  lower  half. 
The  drift  direction  oscillates  wildly  throughout  the  three  magnetically  quite  three-hour  periods 
at  the  beginning  of  the  data  sample  (Kp  is  indicated  below  the  abscissa),  suggesting  the  transit  of 
several  shear  boundaries  through  the  station  zenith.  As  magnetic  activity  increases,  the  convection 
direction  becomes  more  orderly,  with  a  general  direction  of  30°  to  45°  towards  later  times  from  due 
antisunward.  The  drift  velocities  are  low  (150  -  300  msec'l)  during  the  quiet  period,  and  increase 
somewhat  (150  -  600  msec'l)  during  the  disturbed  period. 

Also  of  significance  here  is  the  first  comparison  of  drift  data  with  Interplanetary  Magentic 
Field  (IMF)  data,  shown  in  the  top  half  of  Figure  6.  The  large  fluctuations  in  drift  direction, 
indicating  the  presence  of  shears,  and  the  arc  observations  occur  during  times  when  Bz  >  0,  By  -  0, 
and  Bx  >>  0.  As  Bx  decreases  and  By  increases  with  still  (strongly)  positive  Bz,  the  plasma  drift 
becomes  more  orderly,  but  skewed  by  60°  to  90°  towards  later  times  from  straight  antisunward.  For 
the  steady  8,  <  0  conditons  from  11  to  21  UT,  the  flow  is  steady  antisunward  (skewed  30°  to  later 
times).  As  ftz  becomes  positive  again  after  21  UT,  shears  start  to  dominate  again.  No  clear  effects 
of  By  and  Bx  on  the  drift  direction  (under  Bz  <  0  conditions)  can  be  derived  from  this  short 
measurement. 

Since  for  the  drift  analysis  the  convection  is  assumed  to  be  uniform  across  the  sky  map,  which 
clearly  is  not  the  case  during  the  9  hour  period  of  arc  observation,  the  resulting  drift  directions 
and  velocities  are  not  meaningful  but  merely  give  an  indication  of  the  existence  of  shears.  A  pre¬ 
liminary  attempt  to  separate  skymaps  into  regions  within  and  on  either  side  of  arcs,  using  ASIP 
images  as  a  guide,  shows  evidence  of  shears.  The  detailed  discussion  of  this  analysis  is  beyond  the 
scope  of  this  paper. 


Figure  5.  First  comparison  of  F-regiun  drift  measurements,  taken  on  9/’0  Deccioer  1983  at  Thule,  and 
i nterpl anetary  magnetic  field  (IMF)  data  (for  discussion  see  text).  Also  shown  are  the  Kp  values  for 
the  measurement  period.  y 
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c.  Qaanaaq  -  First  Extended  Drift  Measurements 

Drift  measurements  until  early  1986  had  been  made  by  mannually  selecting  frequencies  and  setting 
sampling  gates,  therefore  requiring  around  the  clock  manning  of  the  Digisonde.  Since  1986  this  pro¬ 
cedure  has  been  automated  for  all  Digisonde  256  stations.  The  automation  is  based  on  the  ARTIST 
trace  identification  software  for  the  Digisonde  (Reinisch  and  Huang,  1982,  Reinisch  et  al ,  1982) 
residing  in  an  IBM  AT  microcomputer,  which  is  part  of  the  Digisonde  256.  The  automatic  drift  program 
also  resides  in  the  ARTIST  computer  (a  report  on  this  important  advance  is  in  preparation).  After 
completing  the  0-trace  identification  the  program  automatical ly  selects  appropriate  frequencies  and 
sets  the  gates.  Throughout  each  measurement  (between  two  ionograms)  the  gates  automatical ly  follow 
the  changes  of  the  echo  height  using  a  two-gate  scheme. 

The  data  (Figure  7)  cover  a  continuous  period  of  almost  6  days  of  drift  measurements,  from 
29  October  to  3  November  1986.  A  cursory  inspection  suggests  at  least  for  the  first  four  days  steady 
antinsunward  drift  direction.  However,  clear  deviations  are  seen  for  the  period  15  -  18  UT  on 
29  October,  and  12  -  16  IJT  on  1  November.  Substantial  direction  oscillations  similar  to  those  dis¬ 
cussed  in  section  b  and  suggesting  the  presence  of  shears/arcs  are  observed,  starting  17  IJT  on 
1  November  under  quiet  conditions,  continuing  more  or  less  unabated  through  the  end  of  the 
measurement. 


Figure  7.  Continuous  F-region  drift  measurements  taken  at  Qaanaaq,  Greenland,  from  28  October  1986 
(21  UT)  to  3  November  1986  (15  UT).  Kp  values  are  shown  between  corresponding  3  hour  tickmarks. 
Superimposed  on  the  presentation  of  the  drift  velocity  is  for  comparison  an  analog  representation 
of  the  kp  time  variation. 

Along  the  time  axis  (abscissa)  we  have  indicated  the  respective  Kp  values.  A  comparison  shows 
that  during  the  more  disturbed  periods  (Kp  >  3)  on  29  October  and  early  30  October,  antisunward 
convection  at  velocities  greater  than  300  msec-1  exists.  This  active  period  is  interrupted  by  a 
3  hour  quiet  period,  (Kp  =  1)  during  which  the  delft  turns  sunward,  in  agreement  with  the  previously 
discussed  satellite  observation  for  Bz  >  0.  A  second  quiet  period  following  a  disturbed  period  with 
steady  antisunward  drift  started  at  12  UT  on  1  November  1986  (Kp  =  1-).  Again  rather  stable  sunward 
drift  is  observed.  From  here  to  the  end  of  the  measurements  the  magnetic  activity  remained  generally 
quiet  and  clear  evidence  of  shears  was  observed  throughout,  even  though  a  short  period  of  moderate 
activity  (Kp  =  3+,  3-)  within  this  quiet  period  tended  to  show  a  somewhat  more  antisunward  organization 
of  the  drift.  The  only  puzz’inq  observations  are  those  taken  on  31  October  1986,  where  Inspite  of 
extended  rather  quiet  magnetic  conditions  the  drift  is  consistently  antisunward. 

The  (horizontal)  drift  velocities  are  shown  above  the  respective  drift  direction  panels.  Super¬ 
imposed  is  the  variation  of  Kp  as  a  function  of  time  in  analog  presentation  (continuous  line).  While 
the  correlation  between  drift  velocity  and  Kp  is  far  from  convincing,  there  is  evidence,  that  higher 
Kn  values  correlate  with  higher  velocites  and  lower  Kp  values  with  lover  velocities  (e.g.  21  IJT  on 
31  October  to-  15  UT  on  1  November).  The  correlation  breaks  down  during  the  extended  period  of  shear 
observation  (2nd  half  of  1  November,  and  all  of  2/3  November  1986). 
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d .  Qaanaaq  -  Drift  Measurements  Under  Sunlit  Conditions 


A  limited  set  of  drift  measurement s  taken  at  Qaanaaq  under  fully  sunlit  conditions  corroborates 
tne  observations  in  darkness.  Drift  measurements  conducted  at  Qaanaaq  on  2y/30/31  July  1986  (data 
not  presented  here'  snow  consistent  antisunward  convection  with  velocities  between  160  to  450  msec"' 
fur  tiie  prevailing  moderately  active  conditions  ( Z  kp  =  26-,  20,  20+,  respectively,  or  Kp  >  2). 

I  jr  measurements  ta*en  under  magnetically  quiet  conditions,  in  full  sunliqht  on  8  April  1986  (Fig¬ 
ure  ,  witn  p,,  -  I-,  i»,  iv,  ?»  for  the  period  of  measurements  (00  to  1?  UT  on  h  April  lyH6  or  22  to 

09  CClT'  tne  convection  is  initially  sunwari,  as  expected  for  quiet  or  (interred)  l>7  >  0  conditions, 
«■>  magnetic  activity  increases  through  Kp  =  1»  to  ?+,  the  drift  direction  turns  antisunward,  skewed 
to  1  a  t  ■  *  r  ti">eS  by  30°  to  6.'°. 

It  is  clear  from  tne  set  of  Qaanaaq  observations  discussed  so  far,  that  Kp  is  at  best  a  coarse 
ordering  .  jrameter  of  tr'e  drift  observations,  especially  since  Rz  >  0  does  not  necessarily  corn— 

S;.  ird  rm  ,1,-t  r  aghet  i  c  conditions ,  as  Figure  o  indicates,  where  a  Jump  from  Bz  =  *3  nT  to  *  9  nT 

v.  amesp  s  to  a  change  in  >  rj  from  Q0  to  30,  contrary  to  expectation  of  continued  quiet  conditions. 
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:',/t  .  i'-regi  >'  drift  observations  from  Qaanaaq,  Greenland,  take  on  7/8  April  1986  taken  under 

■  -i  i -a  ’  I  y  quet  conditions. 

e.  j ,  Labrador  Observations 

As  -ig.,res  1  and  4  indicate.  Goose  Bay  for  the  selected  average  convection  pattern  is  under  tne 
ir*l  , .-rue  )t  the  Si.nward  return  plasma  flow  for  ♦  six  hours  around  CG  midnight,  with  westerly  flow 
expected  prior  to  and.  easterly  flow  past  midnight.  At  other  times  Goose  Bay  is  under  the  influence 
of  o  jrotat  m  g  plasm,!,  measuring  the  lower  drift  velocities  of  the  mid  latitude  ionosphere  (Kent  and 
wr  1  gr  t  .  1  ,  . 


>'  typical  exar.pl e  of  Goose  Bay  drift  measur‘""ents  are  the  data  from  14/15  January  1983  (Fig¬ 
ure  9;,  taken  under  moderately  disturbed  conditions  (Z  Kp  =  16+  and  28+  for  14  and  lb  January, 
respectively).  !>e  convect’O"  in  the  pre-midnight  sector  is  consistent  and  towards  the  west  (sun¬ 
ward;.  Tne  convection  sw’tches  towards  the  east,  moving  smoothly  through  the  expected  southerly 
direction.  A  valid  measurement  indicating  a  quick  switch  back  to  westerly  drift  at  02  AST  inter¬ 
rupts  the  smooth  transition  towards  easterly  flow.  It  should  be  noted  that  the  velocities  are 
considerably  lower  than  those  observed  in  the  polar  cap,  with  premidnight  velocities  measured  at 
250  msec*’1  ,  ana  past  midnight  velocities  at  100  msec"'. 

A  substantial  data  set  (38  full  or  partial  days)  collected  between  1981  and  1985  during  the 
deve ’ opmen t  of  the  drift  technique  is  not  of  the  same  quality  and  consistency  as  the  polar  cap  data 
set.  This  is  due  to  a.)  manual  data  collection,  b.)  interference  with  the  F-layer  measurements  by 
blanketing  Fs  during  magnetically  disturbed  periods,  and  c,)  lack  of  solid  echoes  under  low  foF2 
(trough)  conditions,  coinciding  with  quiet  magnetic  activity  and  small  oval  and  convect ion-pattern 
diameter.  Hut  throughout  the  measurement  set,  consistent  premidnight  westward/postmidniyht  eastward 
drift  direction  is  observed.  Short  15  minutes  to  one  hour  drift  reversals  are  however  not  uncommon. 


especially  in  the  post-midnight  sector.  The  drift  velocities  shown  in  Figure  9  are  typical  for  the 
higher  velocities  observed  in  the  full  data  set,  however  these  higher  velocities  may  be  observed 
after  midnight,  as  well  as  before.  Low  velocities  (<  50  msec'1)  are  more  routinely  observed,  and  are 
mure  typical  for  times  away  from  midnight.  A  thorough  analysis  of  the  data  in  the  context  of  avail¬ 
able  I Mh  data  is  planned,  focusing  specifically  on  the  exact  time  of  midnight  switch,  the  relation  of 
siiortl  ived  drift  reversal  to  IMF  conditions,  and  the  daily  duration  of  measurements  of  distinct  con¬ 
vection  conjitjMS.  This  latter  measurement  may  provide  information  on  the  diameter  of  the  convection 
pa  1 1  ern  . 

r.  Argent  1  a  NAS,  Newfoundland,  Observations 

Our  southernmost  measurements  have  been  taken  at  Argentia  NAS  (57°  CGL).  Here  a  Digisonde  256 
ras  beer,  opr  rat-.  1  to  Support  the  Fast  Coast  OTH-ll  Radar  and  AWS,  unattended  since  summer  of  1985  with 
occasional  visits  hy  a  lucal  engineer  for  data  tape  changes  or  to  reset  the  ARTIST  computer  after  a 
pu«er  outage.  Trie  Digisonde  is  controlled  remotely  via  a  dial-up  telephone  capability,  and  drift 
measurements  in  the  automatic  mode  can  also  be  innitiated  and  supervised  remotely. 

Under  quiet  and  moderate  magnetic  conditions,  Argentia  is  likely  to  remain  under  the  influence 
of  the  correlating  ionosphere,  at  or  south  of  the  plasma  pause  (Fiyure  1).  For  low  magnetic  activity 
i,::  Kr.  <  25),  observed  drift  directions  to  date  do  not  exhibit  any  recognizable  patterns.  Measured 
velocities  are  in  the  50  to  150  msec*1  range.  Figure  10,  however  provides  evidence  that  for  the 
!'  j.jiiet  i  c  a  I  1  y  active  period  of  9  to  12  February  1986  the  diameter  of  the  convection  pattern  expands  to 
b'-irig  Argentia  into  the  area  of  sunward  return  flow.  For  this  station,  magnetic  midnight  occurs  at 


Fur  L-ach  of  tne  3  hour  periods  marked  on  the  time  axes  of  the  three  panels,  the  Kp  values  are  as 

fall u ws : 

9/10  February  panel,  starting  at  18  UT:  5+,  4o,  3+,  2o,  lo,  2o,  lo,  2- 

February  panel,  Starting  at  18  UT:  3o,  3+,  2+,  3-,  2-,  3-,  4o,  5- 

11/  1 7  l.oruary  uariel,  starting  at  18  UT :  4-,  4+,  4+,  3-f,  (ending  06  UT ) 

An  inspection  of  Figure  10  in  the  context  of  the  listed  magnetic  activity  shows  evidence  for 
Sur-njrd  (westward)  drift  at  75  to  150  msec*1  prior  to  local  magnetic  midnight,  with  occasional  velo¬ 
city  increases  to  3uU  msec'1  for  K  >  4.  The  data  on  10/11  February  show  for  less  active  conditions 
[/  <  F ,j  <  4-J  a  switch  between  sunward  (westward)  and  eastward  flow  prior  to  magnetic  midnight,  Sug¬ 
gestive  of  a  sequence  of  contractions  and  expansions  of  the  convection  patttern,  alternating^  shif¬ 
ting  trie  convection  pattern  u r  the  co-rotatiny  plasma  over  the  station.  With  reyard  to  post-midnight 
data,  tne  11  1 ebruary  measurements  show  steady  sunward  (eastward)  flow  with  velocities  between  75  and 
Isj  in  sec'1-  for  enhanced  activity,  while  the  morning  data  on  10  February  show  low  speeds  (  ~  50  m 
and  less  well  organized  drift  direction  under  quiet  condition. 


Figure  9.  F-region  drift  observations  at  Goose  Bay  on  14/15  Janurary  1983  show  the  expected  switch 
of  the  sunward  convection  return  flow  from  westward  to  eastward  near  midnight. 
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Figure  10.  F-region  drift  observations  at  Argentia,  Newfoundland  on  9-12  February  1986  during 
moderately  disturbed  conditions. 


DISCUSSION  OF  DRIFT  MEASUREMENTS 

A  newly  developed  ground  based  technique  to  measure  plasma  drift  at  high  latitudes  has  produced 
data,  which  essentially  agree  with  previous  observations  from  satellites,  and  from  incoherent  scatter 
radars : 

a.  In  the  central  polar  cap  we  find  antisunward  convection  for  Bz  <  0,  or  for  Kp  >  3. 

b.  During  quiet  conditions,  especially  when  following  an  active  period,  sunward  convection 
is  observed  in  the  polar  cap. 

c.  For  extended  quiet  periods  and  Bz  >  0,  shear  conditions  typical  for  the  sunaligned  polar 
cap  arcs  are  observed. 

d.  At  oval/suboval  latitudes  we  observe  sunward  flow  (westerly  before  midnight  and  easterly 
after  midnight),  under  magnetically  disturbed  conditions. 

e.  Argentia  (57°  CGL),  and  to  some  extent  also  Goose  Bay  (65°  CGL)  show  random  switches 
from  westerly  to  easterly  and  back  under  more  moderate  (Kp  <  2)  conditions,  suggesting, 
that  the  station(s)  move  al ternatingly  under  the  influence  of  convecting  or  co-rotating 
pi asmas . 

Important  parameters  of  the  plasma  convection,  namely  the  antisunward  convection  velocity  and 
the  approximate  diameter  of  the  convection  pattern  can  be  estimated  for  at  least  part  of  the  day. 
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DRIFT  MEASUREMENTS  AND  MODELING  OF  THE  HIGH  LATITUDE  IONOSPHERE 


The  importance  of  polar  convection  measurements  lies  in  the  relation  of  the  convection  pattern 
and  the  convection  velocities  (magnitude  and  direction^  to  the  condition  of  the  Interpl anetary 
Magnetic  Field  (l^F),  which  controls  the  polar  cap  potential  tpC .  through  the  solar  wind- 
magnetospheri c  energy  coupling  function  c  (Akasof'i,  1981,  1984,  Reiff  et  al  . ,  1981).  Ionospheric 
modeling  has  shown  that  polar  cap  and  auroral  oval  ionization  and  the  development  of  the  F-layer 
trough  to  the  south  of  the  oval  are  controlled  (location  anil  intensity)  by  the  magnitude  of  the  polar 
cap  potential,  and  therefore  related  to  the  polar  cap  convection  velocities.  We  expect  that  future 
studies  will  lead  to  the  development,  of  ei*her  r  (  IMF)  or  *pC  (IMF)  dependent  high  latitude  iono¬ 
spheric  models  which  can,  in  the  ahsen,o  i'  mpa su remen t s  for  either  parameter,  be  driven  by  polar 
cap  plasma  convection  measurements. 

As  an  example,  we  have  investigated  •  >  “er ts  of  Rz  on  the  electron  densities  in  the  central 

polar  cap.  Using  solar  production  f  nr  *  dependent,  modeling  of  the  ionization  buildup  (and 

decay)  within  a  flux  tube  (for  details  see  Anderson  et  al.,  1987,  this  volume)  as  it  travels  along  a 
trajectory  described  by  a  8.  dependent  rrye^  ion  model  (Heel  is,  1984),  the  time  history  of  the  maxi¬ 
mum  electron  density  (N  maxi  in  flu*  *u‘ies  w'nrh  cross  Thule  at  selected  times  has  been  computed. 

The  effects  due  to  particle  precipitation  have  peon  neglected,  since  the  contribution  to  f*nax  by  the 
cusp  prec  i  pi  tat  i  on  is  small  (  ~2.4  x  in'1  cm-3)  due  to  the  high  cusp  traversing  velocities  of  be¬ 
tween  500  to  1000  msec'l,  (Knudsen  et  al.,  Im//). 

Figure  H  shows  the  trajectories  of  the  flux  tubes  which  cross  over  Thule  at  1400  COLT  (■)  and 
2000  COLT  (A),  in  a  CGL/CGLT  coordinate  system  for  Rz  <  0  and  Bz  «  0  condition,  respectively. 

Other  trajectories  are  also  shown,  hut  are  not  discussed  further.  The  location  of  the  flux  tubes  are 
indicated  in  At  =  30  min  increments;  the  starting  point  of  each  trajectory  is  six  hours  before  the 
time  of  crossing  over  Thule.  The  figure  illustrates  the  large  difference  in  convection  pattern  dia¬ 
meter  and  in  the  convection  velocities  (v  is  proportional  to  the  distance  between  half  hour  marks) 
for  the  two  selected  llz  conditions.  Of  special  importance  is  the  exposure  of  the  relevant  flux  tubes 
to  solar  UV.  The  relation  of  the  flux  tube  trajectories  to  the  terminator  for  December  is  shown  in 
Figure  12.  .  or  Rz  <  0,  the  flux  tubes  arriving  at  Thule  at  1400  and  2000  CGLT,  spend  at  least  4.5 

hours  in  sunlight,  much  of  that  time  at  solar  zenith  angels  of  <  80°.  Due  to  the  high  convection 
velocities,  it  takes  the  flux  tubes  only  between  15  and  30  minutes  to  reach  Thule,  after  crossing  the 
terminator  into  darkness,  a  time  short  compared  to  the  recombination  time  of  the  F-layer  ionization. 

In  contrast,  for  Bz  *  0,  the  respective  flux  tubes  experience  sunlight  for  only  1.5  to  3.5  hours,  and 
only  with  solar  zenith  angles  of  >  85°.  As  a  consequence  of  the  low  convection  velocities,  it  takes 
the  flux  tubes  between  45  minutes  and  1.5  hours  to  move  from  the  terminator  to  Thule.  Figure  13  shows 
the  results  of  the  modeling  of  the  ionization  time  histories,  as  the  flux  tubes  travel  along  the  re¬ 
spective  trajectories.  The  results  indicate  a  difference  of  a  factor  3  (2)  in  N  max  for  ionization 
observed  at  Thule  at  14  (20)  CGLT,  for  the  two  considered  conditions,  Bz  =  0  and  Bz  <  0,  respectively. 


COOIXCTEO  GEOMAGNETIC  L0C4l  time  corrected  geomagnetic  local  time 

CORRECTED  GEOMAGNETIC  LATITUDE  8*<0  CORRECTED  GEOMAGNETIC  LATITUDE 


Figure  11.  Convection  trajectories  for  flux  tubes  crossing  Thule  at  02  CGLT  (  x),  08  CGLT  (O), 
14  CGLT  (□)  and  20  CGLT  (A),  for  Bz  <  0  and  Bz  »  0  conditions,  in  CGL/CGLT  coordindates 
(At  =  30  min).  The  20  CGLT  trajectory  has  not  been  shown  in  the  left  diagram,  for  clarity, 
since  it  follows  closely  the  14  CGLT  trajectory. 
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Figure  12.  The  same  convection  trajectories  as  in  Figure  11,  here  in  geographic  latitude,  local 
solar  time  coordi nates .  Also  shown  is  the  terminator  (  x  =  90°). 


Figure  13.  Calculated  Nnax  values  along  convection  trajectories  crossing  Thule  at  14  and  20  COLT. 


CONCLUSION 

It  has  been  shown,  that  a  chain  of  Digisondes,  using  a  newly  developed  automatic  drift  capabil¬ 
ity,  is  able  to  measure  important  parameters  of  the  polar  convection  patern,  such  as  magnitude  and 
direction  of  the  polar  cap  convection  and  the  approximate  diameter  of  the  convection  pattern.  The 
importance  of  this  new  technique  is  that  it  permits  affordable,  continuous  plasma  drift  measurements. 
With  well  considered  additions  to  the  current  chain,  these  measurements  can  provide  sufficient  infor¬ 
mation  to  lead  to  real  ttme  def Inition/monttortng  of  the  polar  convection  pattern,  or  could  e»  a  mini¬ 
mum  provide  the  time  dependence  and  spatial  extension  for  the  sporadic,  but  large  scale  transpolar 
satellite  drift  measurements. 
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The  dependence  of  the  polar  cap,  oval  and  trough  F-layers  on  the  magnitude  and  structure  of  the 
polar  cap  potential  has  been  well  established  through  modeling  (see  discussion  In  the  Introduction). 
Since  the  drift  velocity  in  the  polar  cap  is  directly  related  to  the  polar  cap  electric  field  (polar 
cap  potential),  this  quantity  is  considered  a  likely  candidate  for  a  driver  of  future  physical  or 
empirical  high  latitude  ionosphere  models.  A  model  study  sketched  out  in  this  paper  shows,  that  for 
different  Rz  conditions  (Bz  <  0  and  Bz  •  0),  the  drastically  different  convection  patterns  and  velo¬ 
cities  (Heelis,  1932)  result  in  substantial  differences  in  the  plasma  densities  in  the  polar  cap. 
This  is  in  good  agreement  with  observations.  We  will  persue  the  outlined  strateyy  in  the  hope  to 
find  ’o  the  polar  cap  plasma  drift  a  new,  measureable  parameter  to  drive  future  high  latitude 
i onosphere  model s  . 
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